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Long distance effects in semi-inclusive B decays 
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We discuss some issues on factorization of long distance effects for semi-inclusive B decay spectra in full QCD 
and in the effective theory. 



1. Introduction 

Let us consider B meson decays into a QCD 
jet plus non QCD partons, such as a real pho- 
ton or a lepton pair. This decay is character- 
ized by three fundamental mass scales, the heavy 
flavor mass rub and the energy Ex and invari- 
ant mass of the jet mx- The threshold region 
is defined as the region where mx <C Ex, that 
is the region where the final observed state has 
large energy and the emission of real soft and 
collinear partons is inhibited. Close to thresh- 
old, in the perturbative calculation, the partial 
cancellation of infrared real and virtual contribu- 
tions produces large logarithms in the ratio of mx 
and the hard scale Q. In the heavy quark effective 
theory (HQET), the heavy quark is replaced by 
a static color charge, leaving unchanged infrared 
exchanges in the heavy meson, and therefore the 
structure of the large infrared logarithms. By in- 
tegrating out the heavy flavor mass in HQET, the 
only remaining scales in the hadronic subprocess 
are mx and Ex and the hard scale Q is identified 
with Ex- Infrared large logarithms at threshold 
occur in the ratio Ex / mx |1I2| . Such logarithms 
have to be resummed in order not to spoil the per- 
turbative expansion. We will elaborate on some 
issues related to the universality of long distance 
effects. 

2. Universality of QCD form factor 

For the semileptonic decay B X^lv, the 
most general differential distribution can be writ- 



ten in a factorized form 
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The distribution is normalized to the radiatively- 
corrected total semileptonic width F. We have 
two short-distance, process-dependent functions, 
the coefficient function C [x, w; a{w mb)] and the 
remainder function d[x,u,w; a^w mb)]. All the 
infrared logarithms are resummed into a univer- 
sal, long-distance dominated, QCD form factor 
a [u; a{wmb)]. The running coupling constant is 
set to the hard scale of the system, Q, where 
Q = w mb = 2Ex ■ 

In the radiative decay B — » Xs 7, a similar for- 
mula holds for the single differential spectrum in 
the hadronic mass 
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= Cr [a{mb)] o [i^; a(mfc)] + 
-I- d_R [t^; Q!(mb)] , 
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where tg — m^ /m'^. The form factor a, being 
the universal form factor, is the same appearing in 
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formula , while the coefficient and the remain- 
der function, dependent on the process, are differ- 
ent. In the radiative decay, it; 1, allowing the 
replacement a{wmh) —>■ a{mh); the distribution 
Q contains a constant coupling a(mb) ~ 0.22. 
Due to kinematics, such replacement cannot be 
done in the semileptonic differential distribution 
The important issue is that the form factor 
a, although being universal, depends on one kine- 
matical variable /m^ in the radiative case, 
while it depends on two kincmatical variables, u 
and w, in the semileptonic case. Such difference 
leads to a natural division into two classes for the 
double and single differential distributions, ob- 
tained integrating the general triple differential 
distribution |2I3I4| . The first class contains 
distributions not integrated over the energy Ex, 
f.i. the single distribution in Ex- they have the 
same form factor as radiative decays. The sec- 
ond class contains distributions integrated over 
the energy Ex, f.i- the single distribution in mx 
or in El : their dependence on the universal form 
factor has been spoiled by the integration, and 
they have the same form factor as radiative de- 
cays only at leading logarithmic order. Since soft 
logarithms signal long distance effects, there is a 
different long distance structure with respect to 
distributions in the first class and to radiative de- 
cay, due to the integration over the hard scale. In 
other terms, while in radiative decays the form 
factor a is fixed at a hard scale Q = mi,, in the 
second class the form factor is integrated on Q 
with a generic weight function (j){Q), dependent 
on kinematics, from low values up to mi,. 

3. Semileptonic spectra 

Let us consider, f.i., the hadron energy spec- 
trum, dT/dEx, belonging to the first class men- 
tioned before. As well known, its shape presents 
the so called Sudakov shoulder, related to the 
occurrence of infrared singularities close to the 
threshold Ex = mb/2. At lowest fixed order 
0{a), the parton process b ulvg, where g 
is a real gluon, contributes to the decay. Above 
the threshold, at Ex = Eu + Eg > mb/2, large 
infrared logarithms, of the form log{Ex ~ mf,/2), 
appear and become singular at Ex mt,/2+. 



Virtual contributions cannot cancel such singu- 
larity, since in the virtual process b ^ ulv the 
energy of the up quark, _E„ = Ex, cannot ex- 
ceed the energy mi,/ 2. The solution is to abandon 
the fixed order calculation and resum into a form 
factor all large infrared logarithms. Resumma- 
tion completely eliminates the singularity, leav- 
ing only an effect in the characteristic Sudakov 
shoulder. 

The single differential distribution in Ex is ob- 
tained by integrating Let us use the adi- 
mensional variables defined in (O and J^J. Since 
there are large logarithms only for w > \, we are 
interested in the resummed formula in that region 
only. We have [2]: 

- C^{a)[l^C2{a) S[w-l;a{mi,)]] 
+ R{w; a), (5) 

where Ci (a) and C2 (a) are two coefficient func- 
tions, and R{'w; a) is a remainder function, van- 
ishing a.t w — 1. They have a perturbative ex- 
pansion in a, whose coefficients are constant in 
the case of the coefficients, depend on w for the 
remainder. The form factor S is the partially 
integrated or cumulative form factor S{u,a), de- 
fined as: 

S{u; a) = I du' <t{u'; a). (6) 

It resums all infrared logarithms, and it is univer- 
sal, being the same as in the radiative decay; at 
any order in perturbation theory the hadron en- 
ergy semileptonic spectra and the radiative one 
have the same infrared logarithms. The form fac- 
tor S{u; a) has an exponential form: 

S = e«, (7) 

with 

cxD n+1 
n=l k=l 

L = log-. (8) 

u 

The situation changes radically if we take a dif- 
ferential spectrum belonging to the second class 
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defined above, f.i. the single differential semilep- 
tonic spectrum in the hadronic mass or the distri- 
bution in the light-cone momentum, normalized 
to mb, p+ = {Ex — \px\)/'mb- After integration, 
such distributions can be rearranged in a factor- 
ized form, with the requirement that the form 
factor contains all the infrared logarithms; how- 
ever, the logarithmic tower will be different with 
respect to the radiative decay |3l4j . 

F.i., in the distribution 1/T dT /dp+, we can 
construct a cumulative form factor Z'p(j)_|_; a) in 
analogy with ©. The related form factor can 
also be exponentiated: 

Up = e^" (9) 
with 

oo n+1 
n=l k=l 

Lp = log—. (10) 

P+ 

The difference in soft logarithms explicitly 
shows up when we compare the numerical coeffi- 
cients Gnk and Gpnk- 

Gp\2 ~ Gi2, 

Gpii = Gil + 0.18, 

Gp23 = G23, 

Gp22 = G22 + 0.089 -0.0046 n/, 

and so on; Uf is the number of active fiavors. 

4. Soft and collinear emissions 

The QCD form factor S can be written as a 
convolution of two functions, each of them corre- 
sponding to distinct sources of large corrections: 
the shape function /, resumming infrared loga- 
rithms due to soft contributions, and a coefficient 
or jet function J . The shape function has the 
characteristics of a quark distribution function 
and it is well defined in HQET. Upon introducing 
a factorization scale /xp we can write 

i:[u; a{Q)] = J{u;Q,ixf) ® fiu^^ip)- (H) 

The coefficient function resums all hard collinear 
logarithms, as well as soft effects between the fac- 
torization scale fjLp and the hard scale of QCD. 



The factorization scale /ip can be identified with 
the ultraviolet cut off of the shape function in the 
HQET. The physical basis of such factorization is 
that soft logarithms are related to longer distance 
effects with respect to hard collinear logarithms. 
In the radiative case Q = m^, and it is a natural 
choice setting — mb as well: 

Simj,) = J{mb) fimb). (12) 

In the semileptonic case, that is not always the 
case. At the end of section (jJl we have rep- 
resented schematically the construction of single 
and double differential decay spectra as an inte- 
gration in the hard scale Q weighted by a generic 
function (j){Q), dependent on kinematics. That 
simplification is helpful here as well. Let us write 

/ dQ^{Q)a{Q) = 
Jo 

= / dQc^iQ) J{Q,fiF) fifip). (13) 
Jo 

In order to avoid substantial soft effects in the jet 
factor, one has to take /ip of order Q. Since Q 
is integrated over, also /ip has to be changed and 
one has to know the shape function as a func- 
tion of /ip. This is to be contrasted with the 
radiative case, in which it is sufficient to know 
the shape function at the single point /ip = to;,. 
Let us stress that it is essential to factorize all 
the soft logarithms in the shape function. One 
can then replace the perturbatively evaluation of 
the shape function with a non perturbative one 
(f.i. lattice QCD), including also non perturba- 
tive soft effects, such as Fermi motion of the b 
quark inside the B meson. 

5. Long distance effects 

Quite often, non perturbative long distance ef- 
fects are included in a phenomenological way by 
convoluting the perturbatively calculated spec- 
trum with a non perturbative structure function. 
Being non perturbative, the latter is generally not 
computed, but, f.i., parameterized and fitted to 
the experimental spectrum in the radiative de- 
cay and then used into the theoretical prediction 
for the semileptonic decay. This procedure is not 
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without risks, since, as seen before, the correla- 
tion between different decays is not always imme- 
diate, even at the perturbative level. It introduces 
the dominant source of uncertainty in the CKM 
determination of Vut from the semileptonic decay. 

Another possibility is to include such effects 
without convoluting the perturbative spectrum 
with a non perturbative structure function, thus 
remaining within the perturbation theory frame- 
work. Effects due to the motion of the b quark 
inside the meson will be taken into account modi- 
fying the QCD running coupling constant into an 
effective coupling according to a model 0. The 
effective coupling will also help solving regular- 
ization problems present already in the perturba- 
tion theory, that are traditionally cured by means 
of arbitrary prescriptions. They are due to the 
fact that resummation formulas do not exclude 
very low energy kinematical regions, where the 
running coupling constant hits the Landau pole. 
The effective coupling is built in two steps 
First, a dispersion relation is used to analytically 
extend the running coupling constant a to a cou- 
pling a without the Landau pole; a has the same 
physical discontinuities than a and the same high 
energy behavior, but it has a finite limit at zero 
momentum transfer. The second step consists of 
resumming in the effective coupling a the sec- 
ondary emissions off the radiated gluons. The re- 
summation formula is the standard one, but it is 
controlled by the discontinuity of a, instead than 
a: 

= ^ dsDisc^M. (14) 

The effective coupling has an expansion in pow- 
ers of a and its usage can be considered a change 
of scheme. As in the previous sections, we can 
factorize the decay spectra and calculate the uni- 
versal form factor a. The difference is that now 
in the perturbative calculation the effective cou- 
pling is used. Factorization is easily performed in 
the Mellin space: 

crjv(a) = [ {l-u f-^a{u; a)du. (15) 
Jo 

The form factor has an exponential form in iV- 



space: 

aNia) = e^"("\ (16) 
where the exponent of the form factor reads: 

JO y 

/ [a(fci)] + B [aiQ'v)] 

+ ^[a(gV)]}- (17) 

This formula is the same as the standard resum- 
mation formula, where effective couplings and 
functions have replaced the standard ones. The 
fmictions A{a), B{a) and D{a) have expansions 
in powers of the effective coupling, and are ob- 
tained by matching order for order in a the stan- 
dard formula with (|17|l . 

The form factor in momentum space is ob- 
tained by inverse transform: 

a{t-a)=(^ ^(l-t)-^ aNia), (18) 

where the constant C is chosen so that the inte- 
gration contour in the iV-plane lies to the right of 
all the singularities of fTAr(a). As anticipated, no 
prescription is needed because (Tjv(Q!) is analytic 
for ReA^ > 0. 

6. Comparison with data 

In Fig. ^ the invariant hadron mass distribu- 
tion for the radiative decay, dVr/dmx is com- 
pared with experimental data from the BaBar 
collaboration 6 . The theoretical curve, calcu- 
lated in the model described in the previous sec- 
tion 5 , shows a good agreement with data. The 
model has no free parameters; however, one can 
choose to best fit the theoretical curve to data 
by using as free parameter a{mz), and obtain 
ctsiniz) = 0.123 ± 0.003. Let us observe that 
the data with mx < 1-1 GeV are not represen- 
tative; they show the K* peak, which cannot be 
reproduced in this approach. 

Another good agreement is obtained by com- 
paring the photon energy spectrum with data 
from the Cleo Collaboration 7 , as shown in 
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Figure 1. Comparison between theory and Babar 
data for tlie B — > X^^i invariant hadron mass distri- 
bution at aijnz) ~ 0.123. 



Fig. m In the B rest-frame, t = rnj^/m^ = 
1 — 2E-y/mB- The photon energies are however 
measured in the T(4S') rest frame, in which the B 
mesons have a smah, non-relativistic motion. In 
order to model the Doppler effect, it is possible 
to convolute the theoretical curve for — com- 
puted with a B meson at rest — with a normal 
distribution of = 100 or 150 MeV. 

Among semileptonic decays, let us consider the 
invariant hadron mass distribution, as shown in 
Fig. |31 in comparison with data from the BaBar 
collaboration Points with mx < 400 MeV, 
which are dominated by the tt peak, are dis- 
carded, as well as points with mx > 2.6 GeV, 
which give basically no information on the signal. 
Once again, there is a good agreement between 
theory and data. 

The only distribution where the agreement 
with data is less satisfying is the electron spec- 
trum in semi-leptonic decays, as shown in Fig. 0] 
where data from Belle collaboration have been 
used |5] . Theory predicts a harder spectrum, with 
a broad maximum around 2.1 GeV, while data 
peak at lower energies. More accurate data, and 
possibly perturbative calculations in the complete 
NNLO approximation, are needed to clarify such 
discrepancy. 



CLEG 
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Figure 2. Comparison between theory and Cleo data 
for B Xsj photon spectrum. Dotted line (blue): 
asimz) = 0.118 and a-y = 100 MeV to model 
the Doppler effect (see text); continuous line (black): 
asimz) = 0.117 and = 150 MeV. 
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Figure 3. Comparison between theory and Babar 
data for invariant hadron mass distribution in 
semileptonic decays for as{mz) = 0.119. 
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Figure 4. Comparison between theory and Belle data 
for the electron spectrum in semileptonic decay for 
as{mz) = 0.135. The data and the theory are nor- 
malized to one in the charm background free region 
< < 0.125. 
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